Cooperative Lattice Dynamics and Anomalous Fluctuations of Microtubules 
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Microtubules have been in biophysical focus for several decades. Yet the confusing and mutu- 
ally contradicting results regarding their elasticity and fluctuations have shed some doubts on their 
present understanding. In this paper we expose the empirical evidence for the existence of discrete 
GDP-tubulin fluctuations between a curved and a straight configuration at room temperature as 
well as for conformational tubulin cooperativity. Guided by a number of experimental findings, we 
build the case for a novel microtubule model, with the principal result that microtubules can spon- 
taneously form micron size cooperative helical states with unique elastic and dynamic features. The 
polymorphic dynamics of the microtubule lattice resulting from the tubulin bistability quantitatively 
explains several experimental puzzles including anomalous scaling of dynamic fluctuations of grafted 
microtubules, their apparent length-stiffness relation and their remarkably curved-helical appearance 
in general. We point out that tubulin dimers's multistability and its cooperative switching could 
participate in important cellular processes, and could in particular lead to efficient mechanochemical 
signalling along single microtubules. 

PACS numbers: 87.16.Ka, 82.35.Pq, 87.15.-v 



INTRODUCTION 

Microtubules (MTs) are fascinating biological macro- 
molecules that are essential for intracellular trafficking, 
cell division and maintaining the cell shape. They dis- 
play unique elastic and dynamic behavior, inherited by 
their complex self-assembling nanotube structure. Sur- 
prisingly, despite an enormous accumulation of knowl- 
edge about their structure, the MTs static and dynamic 
properties have challenged all attempts of a fully coher- 
ent interpretation. In this paper we develop the line of 
thoughts leading towards a new theory that provides a 
more comprehensive understanding of these MT's prop- 
erties. The fundamental result is that stabilized MTs 
spontaneously form large scale superhelices of micron size 
pitches and diameters. The MT's super-helical structure 
turns out to be a consequence of a cooperative inter- 
action between its individual subunits that can sustain 
several stable curved conformations. Cooperativity of 
fluctuating internal degrees of freedom in combination 
with the cylindrical MT symmetry lead to a helical state 
with very unique characteristics in the world of macro- 
molecules: MTs are helices that are permanently but co- 
herently reshaping -i.e. changing their reference ground 
state configuration- by thermal fluctuations. In partic- 
ular when clamped by one end MTs undergo an unex- 
pected zero energy motion. As we will see, this could 
be the key for a consistent interpretation of certain chal- 
lenging experimental results not captured by the conven- 
tional scenarios and models. It is worth remarking, that 
the large majority of experiments probing the above men- 
tioned properties are in-vitro experiments on MTs with 
stabilizing drugs. We stress this point, as the experi- 
mentally prevalent presence of stabilizing agents is not 
innocent and could modify MT's properties. Neverthe- 



less there are reasons to believe that the theory developed 
here might be valid also for non treated MTs. 

Before embarking on the road to a "polymorphic MT 
theory" , we first review the conventional understanding 
of the common MT properties as well as some key exper- 
iments which will be our necessary guides towards the 
model we propose. This paper which extends and deep- 
ens a previous short presentation of the idea of polymor- 
phic MTs [1] is written in two parts. In the first part, 
we provide conceptual and graphical explanations of the 
ideas behind the model and investigate consequences of 
the here developed polymorphic MT model. We hope 
that this part is self consistent and didactic enough, so 
that a general reader can grasp the basic underlying 
ideas. In the second more technical part, the mathe- 
matical model is developed and quantitative results are 
presented. The details and derivations are left to an ex- 
tensive appendix. 



SHORT REVIEW OF WHAT WE UNDERSTAND 
ABOUT MICROTUBULES 



Microtubules are cytoskeletal protein filaments of eu- 
karyotic cells fulfilling different structural and mechan- 
ical functions in the cell: MTs act as "cellular bones" 
strongly influencing the cell shape, constitute the main 
routes for molecular motor mediated intracellular cargo 
transport jH, 0] and perform other important tasks like 
stirring the cytoplasm | J . Besides, MTs play a central 
role in the assembly of the mitotic spindle during cell di- 
vision and are at the heart of the functioning of cilia and 
flagella H . This versatility of MTs in a variety of bio- 
logical functions mainly relies on their unique high stiff- 
ness and on their dynamics of assembly and disassembly. 
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The high rigidity of MTs (similar to hard plastic) is due 
to their structure that is known in exquisite detail from 
3D electron microscopy reconstructions 0, 0] : MTs are 
hollow tubes whose walls are formed by assembly of a 
variable number N of parallel protofilaments (PFs) . The 
PFs themselves are built by head-to-tail self association 
of the a/3-tubulin hcterodimer protein subunit (yielding 
a polarity to MTs) whose structure has been resolved by 
electron crystallography [1, 0] . 

In vivo, MTs most commonly appear with 13 PFs 
[Io| (although there are exceptions depending of the cell 
type), whereas in vitro a variety of structures with PF 
number ranging from N — 9 to 16 were observed [Tll.ll2^. 
Transitions of different lattice types (mostly with the 
gain or loss of one PF) within a single MT are also fre- 



quent [13H15I ]. The MT lattice can accommodate for 



these different structures by twisting the PFs around 
the central MT axis although this process is energetically 
costly. The typical lattice twist repeat lengths (pitches) 
are : +3.4, +25, -6/zm for N = 12, 13 and 14 PF MTs 
resp ectively, with +/- denoting right/lefthanded twist 
[lli - lla ]. As we will see the lattice twist is an essential 
property of the model that we will propose later on. Any 
deviation from the most frequent, energetically favorable 
and thus less twisted configuration N = 13 implies an 
internal prestrain in the MT lattice. The latter stress 
can locally deform the end portions of the lattice or even 
destabilize it [l6j |. 

Another internal prestrain in the MT lattice is believed 
to be caused by the tubulin subunit which upon incorpo- 
ration into the lattice hydrolyzes a bound GTP molecule 
converting it quickly into a GDP-tubulin form which has 
the tendency to form a curved state - curling radially 
away from the axis 17, [l8j]. The constraint imposed 



by the MT lattice however maintains the GDP-tubulin 
dimers in a straight unfavorable state - in turn trapping 
mechanical prestrain. In this conventional view MTs are 
seen as internally prestrained but intrinsically straight 
Euler beams. The GDP-tubulin prestrain is also believed 
to trigger rapid depolymerization called the polymeriza- 
tion "catastrophe" [l9(. The stability of MT is regu- 
lated either by the presence of a thin layer of yet un- 
hydrolyzed GTP tubulin dimers at the growing MT end 
(the so called GTP-cap-model, [13, E0] ) or by the bind- 
ing of MT-associated proteins (MAP) or of drugs such as 
taxol. 



WHAT WE DON'T UNDERSTAND ABOUT 
MICROTUBULES 

Mechanical properties of stabilized MTs 

To investigate the mechanical properties of the MT lat- 
tice the presence of the polymerization dynamics is often 
an obstacle. As previously mentioned it can be switched 



off by stabilizing agents like taxol or MAPs. In the bulk 
of the available in vitro studies taxol stabilization has 
been the experimental method of choice for investigat- 
ing the elastic properties of MT. It is believed that taxol 
maintains tubulin dimers in an approximately straight 
conformation and thus prevents depolymerization j22| — 
[24]]. However, direct EM investigations of single taxol- 
GDP-PFs 25| reveals a more complex and interesting 
picture. A taxol stabilized single PF can in fact coexist in 
several conformational states with comparable free ener- 
gies: a straight state Kp F « and a weakly curved state 
with intrinsic curvature K^p « l/250nm (see Fig. la). A 
third highly curved state with k p c f « l/20nm addition- 
ally appears after longer observation periods. Notably 
Elie-Caille et al. [25| pointed out a cooperative nature of 
the straight to curved transition within single PFs. These 
important findings -several conformational tubulin dimcr 
states and cooperative interaction between them- will be 
one central ingredient of our model later on. 

Going from a single protofilament to the whole tube, 
a central mechanical property of interest often measured 
for stabilized MTs, is its persistence length defined as 
l p = B/ksT. Here B oc Y stands for the flexural rigidity 
which is for an isotropic beam proportional to its elas- 
tic Young modulus Y. The persistence length is the 
length scale characterizing the filament's resistance to 
thermally induced bending moments. Several experimen- 
tal approaches have been developed to measure bend- 
ing stiffness of taxol stabilized MTs. One method con- 
sists in measuring MT's thermal shape fluctuations via 
dark-field microscopy [2f| |2tJ or fluorescence light mi- 
croscopy (2^]. Alternatively in several other experiments, 
l p has been determined by applying controlled bending 
forces via hydrodynamics 27 1, optical tweezers [29h32| 
and atomic force microscope tips 33j. Interestingly the 
authors in [26| observed that stabilized MTs are not per- 
fectly straight, but contain some form of quenched cur- 
vature disorder which needs to be subtracted from the 
measurements. In the same vein, the observations from 
dark-field images of the thermal fluctuations of the free 
end of axoneme-bound MTs, show that taxol stabilized 
MTs adopt a three-dimensional helicoid structure [53] 
(we will return to this point later on). 

The persistence length obtained from these different 
experiments is set in the range of 1 to 8 mm with some 
of newer studies going down to 0.1mm (cf below). For a 
standard semiflexible polymer (like e.g. actin or DNA) 
we expect l p to be a material constant, in particular in- 
dependent of the filament's length. However, for MTs, 
the experimental l p data are not only highly scattered 
but extremely confusing on this point. That the persis- 
tence length could indeed depend on the MT length was 
first mentioned and observed in references 31,K12| and 



confirmed by other experimental measurements - probing 
either the thermal movement 34| or the active bending 
deformations by electrical fields of individual taxol stabi- 



lized MTs gliding on a kinesin-coated surface 35j, 136J . In 
particular these techniques gave for short MT segments 
with submicron length a persistence length between 0.08 
and 0.24 mm. 

This intriguing "length dependent stiffness" was also 
investigated by Pampaloni et al who measured the lateral 
fluctuations of MTs grafted to a substrate 37j. These 



authors found a l p falling within a range of 0.11 to 5.04 
mm for MT lengths varying from 2.6 to 47.5 ^m - with 
a strong linear length-Z p correlation. A similar experi- 
ment done in [38| measured the longest relaxation time 
Tmax (L) of MTs for various MT lengths L. It was found 
that MTs exhibit unusually slow thermal dynamics com- 
patible with r max oc L 3 (cf. Fig 9) - in sharp contrast to 
standard semiflexible filaments with r max oc i 4 . In [39j |. 
the relaxation time for long MTs (L > 10/im) extracted 
from a 2-D shape analysis of taxol stabilized fluorescent 
MTs shows an anomalous dynamics on short scales as 
well. Furthermore it has also been experimentally found 
that the rigidity of MTs depends on their growth velocity 
[4Q3. All these experiments measuring the bending stiff- 
ness and the bending dynamics brought the community 
to the conclusion that the " beam of Life" cannot be seen 
as a simple Euler beam. Its complex internal structure 
should determine its elastic and dynamical properties. 
But which internal mode contributes to the now obvious 
elastic complexity, is the important question - a convinc- 
ing answer of which is still missing. 
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Helices and Rings 

Even more intriguing than the issue of MT elasticity 
is perhaps the question : What is the ground state of a 
MT? While the naive answer - a straight rod - would be 
the most accepted view, a number of experiments put 
this mundane simplistic picture in doubt. For instance, 
wavy sinusoidal and circular shapes are frequently ob- 
served when MTs are adsorbed to glass surfaces or con- 
fined between them. In this confined case the Fourier 
mode analysis of MT deformations systematically reveals 
that a few discrete modes have a larger amplitude than 
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the fluctuations around them [26j, [39J, |40j , cf . also 
-supplementary material. This is a strong hint towards 
the presence of some type of " frozen in" curvature - dy- 
namically quenched on experimental timescales. Like- 
wise in motility assays, it is often seen that when MTs 
glide over a surface coated with molecular motors they 
follow wavy sinusoidal and often circular tracks 42, |43j , 
cf. Figs. lb. In this context, particularly interesting is 
the observation by Amos & Amos 42] of the formation of 
permanently circling MTs (see Figs lb). These unusual 
stable circular gliding structures persisted for many cy- 
cles and occasionally straightened again. As written by 
the authors [42j], this suggests that "an intact tubular 
polymer is capable of holding more than one conforma- 



FIG. 1: The empirical evidence for tubulin bistability: (a) A 
single taxol stabilized protofilament can coexist in a straight 
kpf ~ and a slightly curved Kpp ~ l/250nm state (repro- 
duction from [25|]) (b) Taxol stabilized microtubuls in glid- 
ing assay experiments can switch to a stable circular state 
and move on circular tracks (from [43|[43|). Microtubules are 
occasionally observed to switch back and forth between the 
circular and straight states, (c) Taxol stabilized microtubules 
form a three-dimensional helicoid structure with a Vb\xm pitch 
(from [27 



tional state without the help of an external force" . From 
EM images it was inferred that the underlying mecha- 
nism stems from the balance of individual taxol-GDP- 
tubulin dimers between two or more different stable con- 
formations [441 ] . 

An even more clear hint towards the real nature of 
the "frozen in " curvature was -as already mentioned- 
discovered by Venier et al 27| (Fig.lc) who described 
stabilized MTs as wavy periodic shapes with a half period 
of about 7 — 8fim. This observation combined with the 
fact that MTs often went out of focus, lead the authors 
to " suggest that taxol-treated microtubules may adopt a 
three-dimensional helicoid structure of 15/im pitch". 

Therefore it seems that the MT curvature is a persis- 
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tent attribute attached to its lattice. Any serious attempt 
to fully understand the complexity of MT elasticity can 
not avoid the question of its ground state. Indeed under- 
standing fluctuations around a particular state will be 
futile as long as the origin of the state itself remains in 
the dark. 



lar hinge at the MT clamping point. On the other hand 
artifacts that could trivially lead to a "hinged behavior" 
(like loose MT attachment and punctual MT damage) 
were specifically excluded in experiments 37, 38j |. 

Facing all these obstacles it becomes increasingly clear 
that the solution to all puzzles requires a cut and a rad- 
ically different hypothesis. 



The Soft Shear Model 

A first theoretical attempt to cope with some as- 
pects of the described MT mechanical complexity was 
the "soft shear model" (also called "Timoshenko beam 
model", or "anisotropic composite material model") 



[3J, |31 |45j, |46(. In this model the MT is consid- 
ered as an anisotropic fiber-reinforced material [3^, [33] 
with the tubulin protofilaments acting as strong fibers 
weakly linked with easily shearable inter-protofilament 
bonds. Some specific equilibrium statistical and mechani- 
cal properties of that model were investigated in 45 , 4(| . 
An interesting peculiarity and inherent consequence of 
this model is that any local lattice deformation gives rise 
to a long distance curvature relaxation [46[ and can lead 
to a long range interaction along the MT contour. This 
aspect of the "soft shear model" (SSM) is in phenomeno- 
logical agreement with cooperative deformations induced 
by enzymes like katanin. Furthermore this model pre- 
dicts a length-dependent persistence length which ap- 
proximately resembles the measured behavior 33, |38| . 
However in detail it suffers a number of difficulties and 
inconsistencies in particular : 

-The ground state of SSM is straight - in conflict with 
the helical ground state observation (27| . 

- The SSM does not allow for lattice multistability as 
observed by Amos& Amos [Hj]. 

- The predicted value of the shear modulus is extremely 
small [37|,|3l, 51 (10 5 - 10 6 times smaller than the Young 
modulus) . This would imply very strong shearing in bent 
microtubule structures. This however is unsupported 
by other experimental evidence. Indeed observations of 
straight and highly bent MTs shows that bending does 
not significantly modify the relative position of the inter- 
protofilament bonds [471 ] . 

-The dynamics of clamped MTs [33, |38| does not come 
naturally out from the SSM. To reach agreement and 
fit the experimental dynamics the shear model needs to 
introduce an ad hoc internal dissipation of unclear origin. 

- For short MTs (< 4/xm) that model is very far off 
and disagrees with the "plateau" region of the bending 
stiffness vs length relations ( cf. Fig 3 in [Hj]). 

A careful reanalysis of clamped MT experiments, cf. 
Figs. 2,3 in 33, [38|, reveals two features not captured 
by the SSM : the persistence length scales for large 
L approximately as ~ L (without signs of saturation) 
while the relaxation times scale as ~ L 3 . This exotic 
behavior naively suggests the presence of a limited angu- 



IDEA OF POLYMORPHIC MICROTUBULE 
DYNAMICS 

The new scenario proposed here is based on the hy- 
pothesis of cooperative internal MT lattice dynamics. 
The two central assumptions of our model are as follows: 

(I) The taxol-GDP-tubulin dimer is a conformationally 
multistablc entity and fluctuates between at least 2 states 
on experimental time scales : a straight and an outwards 
curved state (Figs, la and 2a). 

(II) There is a nearest-neighbor cooperative interaction 
of tubulin states along the PF axis. 

Note that assumption I is very different from the con- 
ventional picture where GDP-tubulin has only one ener- 
getically favorable (curved) state. We will show that a 
model based on I and II straightforwardly leads us to the 
very origin of MT (super) helicity and provides a coher- 
ent explanation for static and dynamic measurements in 
thermal fluctuation experiments. 

In contrast to the soft shear model, the present model 
is elastically isotropic but the monomer curvature is 
bistable. As we will see, the ground state in this new 
model is a highly degenerate 3 dimensional helix fluctu- 
ating between many equivalent configurations. 



Conformational Symmetry Breaking and Helix 
Formation 



In this chapter we will give a simple pictorial panorama 
over the consequences of assumptions I and II. What 
happens when tubulin dimers obeying assumption I are 
trapped in the circularly symmetric MT lattice? Start- 
ing from assumption I we imagine that the straight and 
curved GDP-tubulin state have a certain energy differ- 
ence AG > 0, with the curved state being slightly more 
favorable. The strict preference for the curved state is 
however only true if the tubulin dimers are free i.e. not 
confined to the lattice. The situation becomes more inter- 
esting once they are incorporated into the lattice. Obvi- 
ously the outwards bending tendency of the curved state 
is in conflict with the geometry of the lattice. Switching a 
dimer on one MT side will frustrate the dimers on the op- 
posite side of the tube and prevent them from switching 
at the same time. On the other hand the direct neighbors 
of the curved dimer (on the same MT side) will profit and 
switch easier to the curved state, as the lattice is already 
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FIG. 2: Elements of the "polymorphic tube model", (a) The 
GDP-tubulin protofilaments can fluctuate cooperatively be- 
tween two discrete states. The curved a = f state is en- 
ergetically preferred over the straight state a = with an 
energy gain E = —AG. The junction between straight and 
curved states along the same protofilament are penalized by 
a coupling constant E = +J. b) Competition between tubu- 
lin switching energy and elastic lattice strain energy leads to 
spontaneous symmetry breaking: MT bends to a randomly 
chosen direction and assumes a non-zero polymorphic order 
parameter P. The energy becomes invariant with respect to 
an arbitrary rotation of the polymorphic phase angle (j>. 



slightly "pre-bent" in the correct direction. This peculiar 
interplay of negative and positive interaction gives rise to 
a clustering of curved dimer states into a single block on 
one side of the tube cf. Fig. 2b. 

This immediately raises the question about the orien- 
tation of this curved dimer block. Which MT side will be 
selected, and in which direction will the MT overall bend, 
can only be decided by the process of spontaneous sym- 
metry breaking. That is, if the ground state is a curved 
dimer block, it will be a highly degenerate state (see Fig. 
3) . In turn it can be expected to thermally move through 
the lattice at next to no energy cost (apart from some 
friction) ! This is the most essential feature of the present 
model. 

So far we have considered only a single MT cross- 
section. If the assumption II (cooperativity) would not 
hold the curved state blocks would pick their sides at 
each section completely independently. The tube would 



locally curve in random uncorrelated directions and the 
effect of tubulin multistability would stay essentially in- 
visible on the larger scale. However according to assump- 
tion II the blocks become correlated in orientation and 
prefer to stack on the top of each other. This then leads 
to macroscopically a bent - in fact a circular MT - if the 
PF were not twisted around the central MT axis (see 
Fig. 2b). This provision brings us to the last interesting 
point. As already mentioned MT lattices are generically 
twisted i.e. their PFs are not strictly parallel. With a 
cooperative interaction along the PFs - which are now 
twisting around the tube axis - the final product of as- 
sumption I and II will be a long pitched helix! The pitch 
of the helix should coincide with the lattice twist repeat 
length : +3.4, +25, ~6^m for N = 12, 13 and 14 PF MTs 
respectively. The created " polymorphic helix" , however 
will not be unique and will be able to reshape between 
its N indistinguishable orientation states. 




FIG. 3: The straight state of the microtubule becomes un- 
stable and forms a spontaneous bend with fixed curvature 
pointing towards a randomly chosen direction. The micro- 
tubule can assume one of the N degenerate ground states and 
switches between them at no energy cost - the effective po- 
tential has a shape reminiscent of a Mexican hat. 

When we graft one end of such a poly morphic helix 
onto a surface (as e.g. performed in 37|,|38() the helix will 
still be able to switch between the equivalent orientations 
and perform a motion that we call "wobbling" (see Fig. 
4) . It is exactly this type of motion that can give rise to 
the static and dynamic effects measured in 37|, 



To see this we can approximate the movement of a 
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clamped polymorphic helix that is switching between its 
equivalent ground states with a "rigid conical rotor" (see 
Fig. 4). For such a rotor the transverse displacement p 
of the MT end grows linearly with its length L. Using 
the naive definition of persistence length l p « L 3 /3 (p 2 ) 
(where (} is the ensemble average) and the fact that 
(p 2 ) cx L 2 this apparent persistence length becomes pro- 
portional to the length l p oc L. This scaling (cf. Fig 11) 
is in agreement with the experimental results in [371 |38( 
giving us the first hint that the model is on the right 
track. Appendix E comments on the robustness of this 
conical hinge-like motion against a limited local hin- 
drance of the wobbling mode due to the adsorbed part 
of the grafted MT. 

A further encouragement comes from the study of the 
dynamics of the model. Making again the approxima- 
tion of a rigid conical rotation (induced by wobbling) the 
observed unusual scaling of the longest relaxation time 
r max cx L 3 can also be easily understood. In fact a slen- 
der object of length L rotated along a conical surface has 
a friction constant £ rot oc L 3 . In turn the longest relax- 
ation time is given by diffusional equilibration time on 
the cone, i.e. r max cx £,rot/kBT cx L 3 (cf. Fig. 12). It 
turns out that the model correctly predicts both the ex- 
ponent and the prefactor of the experimental relaxation 
time. 




FIG. 4: Clamped polymorphic microtubule with intrinsic lat- 
tice twist attached at its end to a substrate 37^ •'!*) performs 
a peculiar movement. It forms a polymorphic helix with N 
degenerate ground states and switches between them at no 
energy cost. The approximately conical motion with opening 
angle a leads to anomalous lateral fluctuations (p 2 ) ~ a 2 L 2 
- radically different from all other semiflexible filaments. 



In addition it is further reassuring that the helical 
ground state(s) of the model provides a rationale for the 
observation of MT helices by Venier et al. 27 1 . 



In Vivo Implications of Microtubule Polymorphism 

The bulk of observations cited here is made in vitro 
on taxol stabilized MTs. It is known that taxol inhibits 
MT disassembly by maintaining the tubulin dimer in a 
state that strongly favors polymerization. But from the 
experiments reported above it is also clear that taxol 
does not suppress all tubulin conformational changes : 
taxol-GDP-tubulin dimer has multiple stable conforma- 
25j. How do these in vitro findings relate to MTs 



tions 



It is a common empirical observation in many in vivo 
systems that despite their high bending stiffness MTs are 



often seen curved or highly wavy on micron scales [39, 48 



|52j. For instance in [48j, highly sinusoidal MTs on the 
periphery of living LLC-PK1 cells were observed. These 
shapes are usually explained as a consequence of mo- 
tor induced buckling opposed by a gel-like environment 
that leads to finite wavelength- buckling [49]. While the 
" buckling in a gel" interpretation is physically appealing 
a closer look at the data in [49[ (in particular the accom- 
panying movie material) reveals an absence of correlation 
in buckling directions of neighboring MTs. This obser- 
vation puts a question mark on the participation of a 
background continuum gel - as in this case the strains 
in the gel would necessarily propagate to neighboring 
microtubules and lead to spatially correlated buckling 
events - in sharp contrast with observations. There- 
fore we are left with a robust phenomenon of sinusoid, 
constant wavelength MTs - however without a definite 
interpretation so far. 

The phenomenology of stable rings 1421 and wavy si- 
nusoid MTs forming in gliding assays [431 ] is strikingly 
similar and visually indistinguishable from the pure in 



vivo observations 



48 



-52j. Indeed, in both situations 
highly curved lattice states of very similar magnitudes 
kmt ~ 1 — 2/irn~ 1 are readily observed. This analogy 
between in vivo and vitro cases suggests that tubulin 
dimers - both in vivo and taxol stabilized (in vitro) posses 
an identical highly curved state (np-p sa l/20nm). This 
highly curved state appears to be activated within the 
lattice only under compressive loads and seems to be a 
universal property of GDP tubulin itself - independent of 
taxol stabilization. 

On the other hand the weakly curved state of taxol- 
GDP-dimer (nf,-p w l/300nm)) is soft enough to be ac- 
tivated by thermal fluctuations. The empirical evidence 
for the weakly curved state in vivo is far less obvious 
than for the highly curved state. The slight deforma- 
tions induced by the former would be more difficult to 
distinguish from other MT deforming effects in living 
cells like motors, polymerization forces, presence of lat- 
tice defects, bundling and microtubule associated protein 
action. However the observed phenomenology of length 
dependent persistence length of MTs growing from cen- 
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trosomes in egg extracts [H3| (no taxol present) resembles 
qualitatively and quantitatively the in vitro observations 



371 l38j . This is one more indication that the dynamic 



MT polymorphism could persist also in vivo. 



At least two possible biological implications of the he- 
licoidal polymorphic MT nature come straight to mind. 
First, a curved or helical beam under compressive load 
responds like a mechanical spring and is therefore much 
softer (in tension, compression and bending) than a 
straight beam. Therefore a network of helical MTs might 
be important for the overall mechanical compliance of the 
cell. A perfectly straight MT buckled by an extracellu- 
lar load would be much more susceptible to mechanical 
failure and depolymerization than a soft compliant helix. 
Second, helical shapes are geometrically (topologically) 
prevented from side-by-side aggregation, can thus evade 
the formation of bundles and instead form lose networks 
with much fewer contacts. It appears that a tuned helic- 
ity (that can be switched on or off) could be a good me- 
chanical control parameter for the formation of different 
cytoskeletal structures. While the bulk of the cytoplas- 
mic MTs is preferably in the lose network state (favored 
by helicity) in occasional situations like the neuronal ax- 
ons straight aligned MTs are required. In such a situation 
the bundling could be triggered by switching the lattice 
to the straight state. Remarkably in the process of ax- 
onal retraction the straight axonal bundle is destabilized 
(and eventually contracted towards the cell soma) by an 
apparent transition of the MTs to a wavy coiled state 
[54j very reminiscent in shape to single wavy MTs from 
cell soma (HEi^. 



In general MT's polymorphism might have other, less 
obvious biological implications that still have to be iden- 
tified. In particular a more speculative possibility is 
that the tubulin's allosteric multistability might also be a 
piece in the puzzle of MT " catastrophes" . A cooperative 
curvature switch might trigger a transition from growth 
to depolymerization. Maybe the most fascinating aspect 
could lie in the possibility of signal transmission along 
single MTs via a long range conformational switch. If 
our model is correct this is the most inherent and dis- 
tinct consequence of the underlying mechanism. 

We conclude here the qualitative description of essen- 
tial ideas behind the polymorphic tube model. In the 
following we switch gears and present in more quantita- 
tive detail the mathematical model. The mathematically 
less inclined reader is invited to fly over some figures, 
comparisons with experiments, maybe pick up additional 
concepts (such as polymorphic defects and their dynam- 
ics) and jump to the perspectives section which will un- 
derline the essential biological consequences. 



THE POLYMORPHIC TUBE MODEL 

Until now the discussion was qualitative and in the fol- 
lowing we build the mathematical model of the polymor- 
phic MT. We will provide quantitative arguments and 
model experimental data, justifying thus a posteriori the 
previous discussion. In this section, we will focus on the 
thermally induced weakly curved state in taxol-stabilized 
MTs and leave the case of mechanically induced highly 
curved MTs for further works. 

We model the GDP-tubulin dimer state by a two state 
variable a n (s) = 0,1 which corresponds to the straight 
and curved state at each lattice site. The n = 1, ...N is 
the circumferential PF index and s G [0, L] is the lon- 
gitudinal position variable along the MT centerline. We 
recall that our model is based on the following assump- 
tions: 

(I) The taxol-GDP-tubulin dimer fluctuates between 2 
states - straight and curved - (Fig. 2a) with an energy 
difference AG > favoring the curved state. The energy 
density resulting from the switching of tubulin dimers (at 
a given MT section) is then given by 



&trans (^) 



-r->_ ,°v>(«) 



(1) 



with b ?s 8nm the dimer length. 

(II) There is an Ising type nearest-neighbor coopera- 
tive interaction of tubulin states along the PF axis with 
an interaction energy J > favoring nearest neighbor 
dimers on the same PF to be in the same state. This 
leads to the interaction energy density: 



J v ^A r 

einter ( s ) = ~T" / , , (2<7„ (s) - 1) (2cr„ (s + b) - I) 
z — ^n— 1 

(2) 

The last term missing in our description is the elastic 
energy density of the MT lattice. For a usual isotropic 
Euler beam the material deformations e are related to the 
centerline curvature vector k via s = —k-t with r the 
radial material vector in the cross-section. For a poly- 
morphic MT, modelled as a continuum material made of 
N PFs (N = 11 — 16), its actual deformation will depend 
on the polymorphism-induced prestrain e po i . In this case 
the elastic energy density of the MT can be written as 



<»>4 



Ri 



271 



(e — Spot) rdrda 



(3) 



where the integration in e e / goes over the annular MT 
cross-section with Ri ks 7.5nm, R a rj 11.5nTO the inner 
and outer MT radii respectively. The prestrain s po i 
is a function of the polymorphic state of the tubulin 
dimers. Its definition requires a decomposition of 
the tubulin dimer in an inner part (facing the tube 
axis) and an outer part (facing from the tube axis 
outwards), cf. Fig 5. We assume that each curved 
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dimer state generates a positive prestrain +Spf on its 
inner part and an equal but negative prestrain —erf 
on its outer part. We can then write £ po i(s,r, a) — 

e PF <y n (s) [l[Ri,R -dp P /2] 0) - I[R -d PF /2,R ] 0)] 

■I[2, n+qoSt 2, {n+1)+qos] (a) where /[.] (x) = 1 if x € [.] and 
otherwise ( Heaviside function ) and dpF is the PF di- 
ameter. The parameter q appearing in the angular part 
of e po i is the natural lattice twist that leads to the proper 
geometric rotation of a PF around the tube axis. This 
parameter is lattice type dependent and takes discrete 
values 27r/g = +3.4^rn, +25/zm, — 6um for N = 12, 13 
and 14 PF MTs respectively [U, [l3j] [ll] . We can 
estimate the prestrain epp from the experimental value 
of the single switched PF's curvature Kpp ~ (250nm) 1 
-measured in Ref.[25[ on single taxol-stabilized PFs - 
to be epp = dppKpp/2 ps 10~ 2 . Collecting all energy 
contributions together the total elastic + polymorphic 
energy of the MT is then given by 



E 



MT 



-) ds. 



(4) 



The ground state within this model is determined by the 
interplay of the first two terms e e i and etrans- The last 
term ei n t er rules over cooperativity and is responsible 
for the suppression of defects in the ideal polymorphic 
order (cf. Fig. 7). A large value of the cooperativity 
constant with J » k^TLfb would imply a defect free 
lattice. However the presence of the latter defects (and 
their motion) are a necessary ingredient for the overall 
rearrangement of the helix as discussed later on. 



the MT cross-section. The other is an external fixed lab- 
oratory frame with base vectors (u x , u y ,u z ). Putting the 
MT along the u z axis direction and considering small 
MT angular deflections we have u z ss u^. In this case 
the two frames are simply related to each other by a 
rotation transformation R(s) given by the internal MT 
lattice twist qo, such that (u x ,Uy) = i?(s)(ui, U2) with 



R(s) 



cos (70s —sm (70s 
sin qo* cos (70s 



(5) 



To rewrite e in a more illuminating fashion, we define two 
important order parameters at each MT cross-section. 
The first one of them is the vectorial polymorphic order 
parameter 



N 



27TO 27T7j\ 

mi cos — — + u 2 sin --- ) a n (s) 



Physically, P - a 2D vector at each local material frame 
section attached to the MT (cf Fig. 2b)- describes the 
asymmetry of distribution - a kind of "polarization" - of 
the dimer states. It acquires a non zero value only in 
the case when the curved and non-curved states are az- 
imuthally separated on opposite MT sides. For instance 
the "all-straight" or "all-curved" PF state correspond 
both to the same value P = 0. Besides the vector P 
we need to define a second (scalar) quantity 

N 



M( S ) = 5> n (s) 



+£ 



PF 
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FIG. 5: Strains and deformations in the polymorphic tube 
model. Each tubulin dimer can fluctuate between a straight 
state (7 = and a curved state a = 1 of intrinsic curvature 
k,pf- The curved tubulin dimer generates a positive prestrain 
+epf on its inner part and an equal but negative prestrain 
—Epf on its outer part with strains related to observed dimer 
curvature via epf = (Ro — RijKpF /2. 

To understand the basic behavior of the ideal helical 
ground state without defects - i.e. PFs are individually 
in a uniform state (either curved or straight)- we ini- 
tially restrict ourselves to the simplified energy density 
e = e e i + etrans- To investigate the MT geometry we 
first introduce two reference frames (Fig. 2b) . One is the 
material frame with base vectors (u±, ^2,^3) attached to 



which counts the total number of dimers in the curved 
state at cross-section s - or in the Ising model terminology 
: the "magnetization". After integration of Eq. [3] over 
the cross-section and some algebra the energy density 
e — e e i + e t rans can be written in a more appealing form: 



B 



(k - k po1 ) 2 + k\ (i 7 Af - sin 2 (7r/Ar)P 2 )] (6) 



with the elastic bending modulus B = -^p (R^ 



with K\ = 



_ (Rg-R.!) 2 



t(ri+r:* 



kfp and a dimensionless parameter 



7 = 



Kpp 2NAG 
bBn\ 



(7) 



For small deviations of the tube axis from the u z direc- 
tion, the polymorphic curvature vector iZ po i in the ex- 
ternal coordinate frame (u x ,u y ) is related to P (in the 
internal frame (t?i,U2)) via the transformation in Eq. [5j 



K P oi = cR(s)P 



(8) 



with c = sin (7r/iV) /ci a geometric proportionality con- 
stant. 
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Phase Diagram 

In the absence of defects the energy expression Eq. [6] 
allows us to determine the conformational ground state 
of a polymorphic MT. To this end, we first resort to one 
further small simplification and make the "single block 
ansatz" , i.e. at each cross-section we assume only a sin- 
gle continuous block of p switched PFs. This ansatz was 
previously used by Calladine and has been proven very 
useful in modelling bacterial flagellin polymorphic states 
55 , [56j • In the ground-state configuration the curvature 
is given by k — K po i(p) whose absolute value obtained 
from Eq. [8] is K po i(p) — Kism^np/N). The optimal 
switched block size p — p* can be determined by mini- 
mizing the second term in Eq. |6] which within this ansatz 
becomes 



Bk\ 



(9) 



This gives rise to an interesting MT phase behavior (cf. 
Fig. 6). The latter only depends on the polymorphic- 
elastic competition parameter 7 from Eq. [7] -that mea- 
sures the ratio between polymorphic energy of tubulin 
switching and the purely elastic cost of this transition. 

For 7 < —1 the chemical switching potential AG 
strongly dominates the elastic energy cost bBn\. There- 
fore switching is highly favorable and all the PFs will be 
found in the state £7 = 1. This gives rise to a straight but 
highly prestrained configuration. 

Analogously for 7 > 1 the bending energy contribution 
is too costly for PFs to switch at all. Therefore in this 
regime the PFs are all in the straight state with a = 
and the MT is consequently straight as well. 

For — 1 < 7 < 1 the situation is more interesting. In 
this interval we have a coexistence of two locally (meta) 
stable MT conformations : the straight tube (prestressed 
or not - depending on the sign of 7) and a curved lattice 
state with p* switched protofilaments. For —7 < 7 < 7 
with 7 s» 0.72 the curved lattice state is the absolute en- 
ergy minimum and the straight state is only metastable. 
Therefore in this regime and in the absence of twist 
(go = 0) the ground state configuration of the whole tube 
would be a simple circular arc section (cf. Fig. 2b). On 
the other hand, the ground state of a microtubule bear- 
ing natural lattice twist qo ^ will be helical (cf. Fig. 
4). 

It is easy to see that a stable helical state as observed 
in [2lJ is only possible for a switching ratio p*/N £ 
[1/4,3/4]. This together with k pf = l/250nm from [HI 
and Kpoi (p*) — K\ sin (wp*/N) \ provides us with a direct 
estimate of the radius of curvature K~ ol pa 9— lAfim. Very 
strikingly the latter range reproduces closely the observed 
MT helical curvatures as estimated from Venier at al. 
work [27 I kZ 1 „„„ .j ~ Wjim - lending strong support to 
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the model. Furthermore, taking the helical state stability 



FIG. 6: The phase diagram of the polymorphic microtubule 
model as function of the polymorphic-elastic interaction pa- 
rameter 7 from Eq. [7J Depending of the magnitude and sign 
of 7 the microtubule can be either in an "all PF switched 
state" (black), "no PF switched state" (grey) or in a mixed 
curved- helical state. Only the latter "mixed state" will dis- 
play net curvature and lead to an observable helical appear- 
ance. 



as an empirical fact (implying that 7I < 0.72) and as- 
suming a typical protein Young modulus Y ps 1 — lOGPa, 
allows us a simple estimate of the transition energy per 
monomer AG ps +1.1 to +11&T — a reasonable range for 
a soft biological object. 

In general, the full energy expression -including the 
cooperativity energy term- Eq. U gives rise to a very 
complex behavior. Here we will focus on some basic new 
phenomena. It turns out that the most remarkable de- 
viation from standard wormlike chain (WLC) behavior 
arises from the fluctuation dynamics of the polymorphic 
order parameter's angular phase that we consider in the 
following. 



Polymorphic phase fluctuations 

Here we introduce a sightly different phenomenologi- 
cal model that simplifies the study of Eq. |4] while still 
reflecting important aspects and physical properties of 
it. In this section we assume as before that the helical 
state is the ground state and consider now the effect of 
the fluctuations around it. To this end we decompose the 
polymorphic order parameter 

P(s)=P (s) [cos (<j> (s)) mi + sin (0 (s)) u 2 ] 

with P (s) being the "polymorphic amplitude" and </>(s) 
the "polymorphic phase variable". The latter one de- 
termines the orientation of the switched block tubulin 
dimers at each cross section with respect to the MT ma- 
terial frame. From Eq. [SJ the centerline curvature with 
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respect to a fixed external (u x ,u y ) frame is then: 

~i? P oi(s) = k [cos (q s + 4> (s)) u x + sin (q s + (j) (s)) u y ] 

(10) 

with ko = cP(s). 

In general both the polymorphic phase <p an d ampli- 
tude P can fluctuate along the MT contour and give con- 
tributions to the polymorphic energy. The phase fluctu- 
ations are induced by creation and motion of polymor- 
phic "double defects", cf. Fig. 7. The double defect - 
a kind of " polymorphic dislocation" - that can be either 
left or right handed - maintains the number of switched 
protofilaments constant while reorienting the direction 
of curvature. At zero temperature the lattice would be 
defect-free, <p would be constant and the polymorphic or- 
der parameter P would strictly follow the lattice twist. 
The phase change </>' = dcj)/ds will deviate from zero if on 
relevant length scales there are enough polymorphic dou- 
ble defects to allow for a reorientation of the polymorphic 
order parameter away from optimum. The double defects 
carry only a limited local energy cost AE = 2 J per de- 
fect and can be easily thermally excited if J < fceT. In 
the approximation of a large number of PFs, assuming 
that (j) can change continuously we can write the phase 
contribution to the energy as 

E P oiU) = ^ £ dscj) 12 (11) 



polymorphic phase stiffness 

= 2J/fc £ 



with the 

k B T^r (2 + e 2J / kBT ) which can be related to the 
density of double defects with energy 2J (cf. Fig. 7 
and Appendix A), giving rise to a new length scale - 
the polymorphic phase coherence length 1$ — C^/ksT. 
For MTs shorter than 1$ we will observe coherent 
helices while on longer length scales the helix softens 
significantly and looses eventually its helical appearance. 

In contrast to the just discussed polymorphic dislo- 
cations which can be easily thermally excited, the vari- 
ation of the polymorphic amplitude P, i.e, change of 
the number of switched PFs is more energetically costly. 
Any deviation of P from its optimum state P* (given 
by the phase diagram) is associated with an energy cost 
E oc {\P\ - \P*\)) 2 ■ I proportional to the length I of the 
region in the unfavorable state, cf. Fig. 7 (see Appendix 
B). Therefore we conclude that on large enough scales 
the polymorphic phase fluctuations will be the dominant 
effect. Based on this and on the observation of stable he- 



lical states [27[ we will in the following assume P = const. 
For small deflections around the z-axis, the unit vector 

tangent to the MT's centerline is approximately given by 
t « (6 X , 6 y , 1) in the laboratory frame (u x , u y ,u z ) where 
6 = (8 X ,9 V ) are the centerline deflection angles in x/y 

direction. The global centerline curvature it — d t /ds 

can then be approximated as it w (0^,0^,0) and the 



Ideal Polymorphic 
Order 




Double Junction 
Defect 




Inglt Defect 



"Random Walk" 
of Double Defects 



FIG. 7: Defects in ideal polymorphic order soften the helical 
states and give rise to "polymorphic dynamics". Single de- 
fects carry a cost proportional to their length, double defects 
only a local energy contribution. The coexistence of left- and 
right handed defects (LH and RH) along the length leads to 
a "random walk" of the polymorphic curvature direction and 
in turn to an effective torsional deformation. 



total MT energy can be written as follows: 
Etot = E pol (0) + E e i (0, 0) • 



(12) 



The first energy term is the polymorphic phase contri- 
bution Eq. [TTJ The second term is the elastic bending 
energy 



E„ 



~it po i) 2 ds. 



(13) 



From Eqs. [TTIfTSI we see that the zero-temperature 
ground state corresponds to <f> — const and 0' = ~it po i 
- that is to a defect-free helix with a pitch given by the 
natural lattice twist qo. At finite temperature, both elas- 
tic and polymorphic fluctuations will be excited so that 
the curvature can be decomposed as 0' — lt po i + 6 e i 

with 6 e i the purely elastic contribution. This gives 
rise to a helical MT shape described by the curvature 

~itpoi + v ei and torsion r ~ <j)' + 1o- The MT lat- 
eral displacements away from the z axis can be writ- 
ten as f(s) = (x(s),y(s)) = f° (6 x (s')),9 y (s'))) ds' that 
for small deflections decouples into elastic and polymor- 
phic displacements ~{?{s) = p po i + p e U where ~J? e i ~ 
Jq ei(s')ds' and ~jt po i ~ po i(s')ds'. The latter can 
also be written from Eq. [10] as 



~tpol{s) = k ds' ds (cos (q s + 4> (s)) it x 

JO JQ 

■4><ft)-ty) (14) 



-sinlgos- 



In Fig. 8, snapshots of configurations of clamped MT 
obtained from Monte Carlo simulations are plotted for 
different concentrations of double defects (i.e. different 
values of 1$) with twist and no twist. It is interesting 
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to remark at this point that based on the symmetry in 
the problem any MT configuration can be rotated around 
the z— axis with no energy cost. This seemingly trivial 
feature - the energetic degeneracy- is in fact the most 
distinctive and unusual property of a polymorphic chain. 
We consider the consequences in the following. 



a 



No Defects (Twisted) 




With Defects (Twisted) 





No Defects (No Twist) 



With Defects (No Twist) 



FIG. 8: The conical "wobbling" motion due to the rotational 
energy degeneracy of the polymorphic MT model : Snapshots 
of Monte-Carlo simulated lattice states, (a) Twisted MTs free 
of defects L/lj, <C 1 and with numerous defects L/lj, > 1. For 
larger number of defects, the helix looses its coherent look, (b) 
For non-twisted MTs the movement has a typical parabolic 
"trumpet-like" shape. 



The Wobbling Mode 

By construction a polymorphic chain as we describe 
it here has a N fold symmetry. Therefore, there are N 
different helical states of different orientations with the 
same energy i.e. N grounds states. This energy degen- 
eracy is also reflected in the continuum model (where 
the N fold symmetry is approximated as continuous) by 
the rotational invariance of E po i (</>) (which depends only 
on cj)'). The broken cylindrical to helical symmetry of 
the straight state is then restored by the presence of a 
" Goldstone mode" (p —> <p + 0o consisting of a rotation 
of P by an arbitrary angle (f>o in the material frame (cf. 
Figs. 2 and 3). This mode comes energetically for free 
and leads to dramatic effects on chain's fluctuations. For 
instance for a MT clamped at one end, this symmetry 
implies that the MT will randomly rotate much like a 
rigid rotor as shown in Fig. 4. Note however that this 
rotation will still be associated with a certain dissipation 
as the system has to go over energy barriers between two 
helical states. This barrier due to the flipping of lat- 
tice states can be overcome at nonzero temperatures by 
the creation of double defects which diffusively propagate 
along the MT and eventually angularly reorient the poly- 
morphic order parameter P. These dynamic phenomena 
and the dissipation mechanisms will be discussed in a 
later section. 



In summary, the zero-energy mode that we will also call 
the "wobbling mode" is an inherent feature of a helically 
polymorphic filament and as we will see now could be the 
culprit causing anomalous fluctuations of clamped MTs. 



Persistence Length Anomalies 

Among several definitions of the persistence length, we 
consider here - for direct comparison with clamped MTs 



experiments 37 , 3|| - the " lateral fluctuation persistence 



length" defined as 



11(8) = (2/3) s 3 /V(s) 



(15) 



with V(s) = (^p(s) J — (p(s)) the variance of p 2 = 

x 2 + y 2 - the transverse displacement at position s and 
the ensemble average. We assume as in experiments 
371 " 



a rigid clamping point at the position s = 
preventing the microtubule from translating and rotating 
at that point. 

For an ideal semiflexible wormlike chain we expect that 
the persistence length I* = lg is a position independent 
and definition invariant quantity equal to the bending 
persistence length Is = B/ksT. (For another more clas- 
sical definition of the persistence length - coming from 
the tangent tangent correlation function, see also Ap- 
pendix C). However for a polymorphic chain the strict 
equivalence of I* and Ib is not correct. To see that, we 
can decompose the polymorphic and elastic fluctuations 
p(s) = Ppoi + el- Inserting this in Eq. 1151 and taking 
into account that for small deflections the two compo- 
nents decouple po i~~p* e i) — leads us to the following 
relation for the persistence length 



(16) 



with the polymorphic persistence length given by l po i — 

(2/3) s 3 /V po i and V po i = (p P oi 2 ) - (p P oi) 2 - The average 
(..) is now performed over the phase (f> governed by the 
energy Eq. [TTJ More precisely, the average of any arbi- 
trary functional A [</)] of the polymorphic phase can be 
performed by first selecting one of the equivalent ground 
states denoted 0o characterized by 0(0) = (f>o and per- 
forming the average over the polymorphic angle distribu- 
tion Eq. [TT1 around the chosen ground state, i.e., 



4>o 



D<t>A 



cxp(-- 



L 



ds(i>' 2 ) 



(17) 

where </> = <f> — <fro (and thus 0(0) = 0) and Z = 

J Z)<j6exp(— y Jq ds}' 2 ) is the partition function. In a 
second step - for a freely rotating polymorphic phase- we 
integrate over the rotational zero mode 0o : (^[0]) = 
J (.A [</>]) \4, d4>Q- This operation correctly takes into 
account the phase fluctuations over (and around) all 
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equivalent ground states related by the transform <f> — s- 
4> + 4>q. The rotational symmetry around the z-axis 
(integration on 4> ) readily implies (p po i (s)) = and 

( x poi) = (lipol) ■ Therefore the polymorphic persistence 
length can be written as 

l pol ( s ) = (l/3) S 3 /(y 2 pol ( S )). (18) 

with 



(y P oi( s )) = J ~2^~ J J { 6 y,voi{si)O yiV oi{s 2 ))\4 10 ds 1 ds2 

(19) 

whose computation (for details cf. Appendix C) leads to 
the following mean square displacement 

(Vvolis) 2 ) = 3(1 + ^2)4 { P l W - e_ ^ P 2 (S) COS (q 0S ) 

-e"^P 3 (s)sin(g s)} (20) 

where Pi(s) are polynomial functions given in Appendix 
C. 

A typical curve of I* vs s is provided in Fig. 9. In 
general it shows three different characteristic regimes de- 
noted I, II and III in the figure : 

I. At short distances to the attachment point s < 
Smin ~ tt/<Zo (half the polymorphic wavelength) the total 
persistence length can be approximately given by 



It w h 



In the limit of very short distances s « lg Kq , the 
polymorphic fluctuations become negligible and are com- 
pletely dominated by purely " classical" semiflexible chain 
fluctuations. Not surprisingly the persistence length co- 
incides then with the classical bending persistence length 
I* (0) = Ib- Starting from Ib, polymorphic fluctuations 
begin to contribute reducing I* that attains a global min- 
imum at s m - m » tt/qq. 

II. For intermediate length values s m i n < s < l^, the 
total persistence length displays a non-monotonic oscil- 
latory behavior around a nearly linearly growing average 



2<?o , 4 9o 



(22) 



This result is worth deeper understanding. A moment of 
thinking reveals that the oscillatory part with wavelength 
2n/qo is related to the helicity of the ground state. At 
the same time the linear growth I* (s) ex a 2 s can be asso- 
ciated with the roughly conical rotation of the clamped 
chain (wobbling mode) which acts as an effective "rota- 
tional hinge" at the attachment point, cf. Fig. 4. The 
sinusoidally modulated rotation cone which builds an ap- 
proximate envelope for the chains' motion has an opening 



angle a which is related to the geometric features of the 
helix a = 2koQq ■ 

III. Finally for very large distances from the attach- 
ment point s 3> 1(f) we expect to recover classical re- 
sults of a semiflexible chain again. Indeed in this asymp- 
totic regime the effective persistence length reaches sat- 
uration with a renormalized constant value I* (oo) = 

l/^+'s 1 ) where 



i 

— K. 

2 



o > 



(23) 



Intuitively the helix looses then its "coherent nature" - 
due to strong variations of <fi' an d elastic fluctuations 
6 e i - and the collective rigid rotational ("conical") mo- 
tion is finally replaced by an uncorrelated segment move- 
ment. Not surprisingly the persistence length becomes 
then length independent again. Curves for different val- 
ues of Irf, are provided in Fig. 10. 



Untwisted MTs 

While there are no completely twist-free MTs and ev- 
ery lattice will have generically a small twist, one can 
still formally study the interesting limiting case qa = 0. 
Note that the large estimated pitch of 13PF MTs is fi- 
nite and in the range of > 25/j,m 11- HI 15 1 while often 
assumed to be approximately infinite. In such an ideal 
case the theory still applies, however the overall behavior 



(21) of I* (L) will substantially change and become much less 
consistent with the linear scaling found in experiments. 
While for q$L > 1 the chain in leading order moves on 
a linear cone (with fixed opening angle a) for q$L << 1 
while still L « 1^ the "wobbling" motion is happening 
on a quadratic cone (a "trumpet shaped" cone, cf. Fig. 
8b) . More precisely for the exceptional case of untwisted 
MTs the polymorphic part of the lateral fluctuations Eq. 

HOI behaves as <y po /(s) 2 ) = 2/3/t 2 ,^ (jPi - e~^P 2 ) with 



Pi (a) = 24Z3_3^s 2 + s 3 , P 2 (s) 



24/' ! 



12l 2 s. There- 



fore for short MTs L « 1$, the lateral fluctuations grow 
with the length in fourth power (y po j(L) 2 ) k 2 L 4 /8, 
whereas for long ones L » 1$, the deviation grows 
cubically, (y po ;(L) 2 ) = 2/3kqZ0L 3 . From this, the per- 
sistence length has consequently two typical regimes. 
For L << 1$ we deduce from Eqs. HH] and [T5] that 

l*(L) » (Z" 1 +3/8/4L)" . The latter expression im- 
plies that long, untwisted MTs appear increasingly softer 
with growing length, and the effective persistence length 
decays inversely I* oc 1/L for L >> 1/(3/8kqZb) and 
reaches a limiting value I* = Ib/ (l + 2Ib^«o) for L >> 
1$, For a MT of L ~ 10 — 20/im we would expect 
I* ~ 10 — 20 fxm - a value 2 orders of magnitude smaller 
than observed in 37, [38| . This decreasing behavior is 



in contrast with observations I* oc L - leading us to the 
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FIG. 9: A typical shape of the effective persistence length 
11 (L) for a clamped microtubule as obtained from Eqs. 1161 
I18l and [20l Most generically the curve displays three different 
regimes: (I) An initial rapidly decreasing regime where poly- 
morphic effects become more effective with growing L (soft- 
ening the chain), (II) a linearly-growing oscillatory regime 
- corresponding to the coherent wobbling movement of the 
clamped microtubule, (III) an asymptotic plateau regime, 
where the helix progressively loses its coherence with grow- 
ing L. In this regime, the behavior tends to that of a classical 
semiflexible chain yet with a renormalized effective persistence 
length given by cf. Eq. 1231 



conclusion that the 0-twist MTs does not constitute a 

and 



significant portion of the experimental data [37 , 
twist is necessarily required for the growth of I* with L. 

Having developed some static consequences of the 
polymorphic MT we now turn to its dynamical aspects. 



Polymorphic Phase Dynamics 

To describe the MT fluctuation dynamics we consider 
the total dissipation functional Pdiss — Pext + Pint which 
is composed of an internal dissipation Pi nt = h^mt J 4> 2 ds 
(with <j) = d<p/dt) coming from the flipping of lat- 
tice states and an external hydrodynamic dissipation 

2 

Pext = ^£-L J p ds associated with the time variation 
of the MT deflection (s,t) = (x(s,t),y(s,t)). We as- 
sume that the friction constant (per unit length) £j_ of 
the helical MT is approximately the friction constant 
£i = Atttj/ (ln(2L/i?) - 1/2) of a long slender body of 
length L, (small) radius R « L moving in a fluid with 
viscosity rj at low Reynolds numbers. The time evolution 
equation of the phase variable 4> (s, t) and the lateral dis- 
placement y (s, t) (and x (s, t)) are given by the coupled 
Langevin equations 



SE t , 



SP dl . 



(24) 



FIG. 10: Comparison of theoretical persistence lengths for 
different values of the polymorphic phase coherence length 
l<t>, (Ib = 25mm, kq = O.OSfim -1 and qo = 0.8/^m _1 ). For 
Icj, = oo, the MT is a (defect free) coherent helix performing 
the "wobbling motion" (as in Figs. 4 and 8a left panel) . The 
plateau regime - where elastic fluctuations become dominant 
over polymorphic- is reached for very long MT only (not seen 
in the Fig.). Finite 1$ reduces the coherent wobbling motion 
- shortens region (II) of Fig. 9 - and the plateau regime is 
reached earlier with decreasing 1$. 



and 



SE t . 



SP d 



5y 



5y 



(25) 



with r^/p the thermal noise terms. In general the lateral 
displacement y(s,t) has contributions from both poly- 
morphic y po i{s, t) w k J* ds' J* sin (<j>(s, t) + q s) ds and 
elastic fluctuations y e i(s,i) « J Q S e i(s' ,t)ds' and the dy- 
namics is highly non-linear. In the regime L >> 1$ where 
the helix looses its coherence one expects to retrieve the 
dynamics of the usual semiflexible filament with r ~ L 4 . 
However in the opposite and physically more interest- 
ing regime L << 1$ a new and different dynamic be- 
havior can be expected. As we learnt from the study 
of the static case the effects of polymorphism become 
more pronounced at shorter lengths. As we have seen in 
this regime, the dominant motion is the wobbling rota- 
tion of a coherent helix on a cone where elastic fluctu- 
ations become negligible compared to polymorphic ones 
i.e. y(s,t) w y po i(s,t). In this regime few polymorphic 
defects L << 1$ are present and the phase can be ap- 
proximated as 4>{s,t) k, cj>o(t) + 6(f> (s, t) .Using this de- 
composition with 8<t>{s,t) « 1 we can expand P ext to 
leading order 



P 



'diss « \l {tint + e«rt) 4>o + O (4 2 ) (26) 



with a external friction constant ^ ex t given by 

(,ext = — ^(2(1 + cos Lgo) -4— + -^hr~ ) ( 27 ) 

<7o Lq 3 
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The evolution of the zero mode <j>o (t) reduces from the 
Langevin equation, Eq. [24j to = — sp s . d j ss + which 

Otpo 

leads to the equation of motion 



dt 



<h (t) 



i 



-IT 



T<t, (s,t) ds 



(28) 



with a friction constant £tot = £,mt+£,ext- Therefore 4>o (t) 
satisfies the simple Langevin equation Eq. [^correspond- 
ing to a simple potential free Brownian motion with its 
mean square displacement given by (see Appendix D for 
a more detailed explanation) 



(W>oW-M0)) 2 ) = ^ 



(29) 



In this limit (wobbling mode dominant) we have a 
roughly rigid helix moving randomly along a cone and 
all the physics is contained in the effective friction coeffi- 
cient £t t (£mt, Ko, <Zcb L) and its dependence on the inter- 
nal dissipation Cint, the helix parameters kq, qo and the 
length L. 

For later comparison with experiments we compute the 
longest relaxation time r given by the auto-correlation 
function (y po i(s,0)y po i(s,t)) oc e~ t/r . Using Eq. [29] a 
short computation (Appendix D) leads to 



(y pol (L,0)y pol (L,t)) = (y 2 pol (L)) e ~*M£) 



with 



VUL)) = § ( 



1— cos(qoL) 

53 



l(<?o£) 



(30) 



the 



static mean square displacement Eq. [20] in the limit 
L « 1$ and with r (L) = L^ to t/k B T the longest re- 
laxation time - proportional to the total friction con- 
stant £i t - inheriting its length dependence in two dif- 
ferent length regimes. For very short lengths L <C l c — 
(3£int£j ) qo K o when the hydrodynamic dissipation 
is entirely dominated by internal dissipation we have a 
linear scaling 



r (L) « L^ int /k B T. 



(31) 



For larger L > l c the wobbling movement through the 
fluid is the dominant source of dissipation and 



T{L)^L^ ext {L)/k B T. 



(32) 



with £ ex t(L) given by Eq. [27J Note that the L depen- 
dence of £ ex t relies also on the L dependence of £_i_(L) 
which for simplicity has been modelled as the friction of 
an ideal slender tube moving in a liquid. A more precise 
(but difficult) determination of £j_ could slightly change 
its variation with L, although not the general trends of 
£eiEt(£). In particular if we assume that £j_ is L indepen- 
dent and that Lqo 3> 1 we have in this regime the scaling 



r(L) 



6 



3k B T 



(«o/9b) L 3 



(33) 



So in summary for r we expect a cross-over from a linear, 
internal dissipation dominated L dependence at short 
lengths to a cubic length dependence given by hydro- 
dynamic friction alone. 



COMPARISON WITH CLAMPED MT 
EXPERIMENTS 



The comparison with experiments [371 138( which mea- 
sure lateral fluctuations of clamped MTs reveals several 
interesting characteristics that are in agreement with pre- 
dictions, cf. Fig. 11. First, the predicted mean linear 
growth of I* (L) agrees with experiments as a single expo- 
nent fit I* ~ L 5 of the data provides 5 = 1.05. Besides the 
linear growth the experimental data reveal a large spread 
of I* data points which seems to grow approximately in 
proportion to the length. This linearly growing experi- 
mental spread is likely linked to the intrinsi c sp read of 
q values of different MT lattice populations 11-13, 1S| . 
MTs with different number of protofilaments will display 
different lattice twists ranging from 50 « 27r/3/irn (for 
12 PF MTs) to <jo < 27r/25^m (for 13 PF MTs). Keep- 
ing in mind the scaling I* (s) cx (JqS (cf. Eq. [25]) we 
would expect more than an order of magnitude variation 
of measured I* while the slope I* (s) / s should display a 
constant spread. 

Second, the data of Taute et al. [38[ (Fig. 11, circle) 
indicate a non-monotonic dependence with systematic 
trends over several consecutive data points. This seems 
phenomenologically well captured by the oscillatory be- 
havior I* (L) from E q. [2"2l On the other hand, the data 
of Pampaloni et al. |37| are definitely much more spread 
and don't allow such a clear conclusion. Therefore the 
non-monotonicity of I* (L) is at present experimentally 
difficult to infer from the two existing experimental data 
sets taken together, however it is consistent with data 
within the error bars. As mentioned the presence of dif- 
ferent lattice populations (even within single MTs) could 
give rise to a large spread of experimental data points 
and an effective " washing out" of the non-monotonic be- 
havior for different lattices within the same statistics. 

Remarkably the experimental data reveal that the 
large length plateau s 1$ where I* would become 
length independent is not reached even for longest MTs 
(~ 50)im). This is in phenomenological agreement with 
the theory - as based on the observed long coherent he- 
lices by Venier et al [27| it would imply a very long 1$. 
The absence of the plateau in Pampaloni and Taute's 
data allows a lower estimate of the coherence length: 
l,p > 55/im which in turn would imply a large coupling 
constant J > Ak B T. The best comparison between the- 
ory and experiments (cf. Fig. 11) gives l B — 25mm cor- 
responding to a rather high Young modulus Y ss 9GPa 
higher than typically reported before 26-32, [3^, |40| . 
However the present value is well within the range for 
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proteins and protein tubes with Y up to 19GPa are re- 
ported in literature [57j . The higher value of the bare 



Young modulus extracted form the present theory should 
also not be a surprise as in previous models all MT con- 
formational fluctuations were interpreted as originating 
from bending deformations alone. In our model - both 
clastic and polymorphic fluctuations contribute with the 
latter being much softer and giving therefore dominant 
contribution. 

The best fitting helix wave length A rj 1 .hii m is clos e to 





11 13 



the expected 6/im corresponding to the twist 
of the 14 PF MT population. This is in agreement with 
the fact that, - in contrast to the in vivo situation- the 
large pitch (A * 25/j,m) 13PF MTs are likely underrep- 



resented in the data [33, |38j . Indeed in vitro studies of 
taxol-copolymerized MTs display a MT population con- 
sisting of a majority of 14 PFs (61%) while 13 PFs (32%) 
are less represented 11- 13, Ell- The in vitro conditions 
therefore strongly shift the PF population away from the 
preferred low twist 13 PF MT towards the highly twisted 
14PF MT. 
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FIG. 11: Effective persistence length Z* as a function of the 
position from the attachment point along the clamped MT 
contour. The experimental data (stars and circles) [37l . [Hj] 
and the theoretical prediction with Ib = 25mm, A = 7.5/im, 
Kg 1 = 18fxm, qolj, S> 1. 

The estimated Ib is larger than in previous studies 
(Ib ~ 1 — 6mm) where however polymorphic fluctuations 
were neglected. This result leads us to the conclusion 
that if polymorphism is partially suppressed one would 
measure much larger effective l p - as in fact observed. For 
example in studies where 2D slab geometry is used (MTs 
between 2 close glass slides) an effective suppression of 
the 3 dimensional polymorphic helices or their reduced 
mobility is expected. A typical observation in such cases 
is an extensive "intrinsic curvature" (of previously un- 



known origin) . Within our theory one could interpret this 
curvature as pinned / quenched polymorphic helices pre- 
vented from free fluctuations by the confinement. These 
effects could in general explain the dramatic variations 
of measured l p values based on the presence / absence of 
polymorphic "softening" in different experimental setups 
and geometries. 

Now, let us consider the clamped MT dynamics as in- 
vestigated in [38|. A careful analysis of Taute et al's 
data reveals a peculiar scaling of the longest relaxation 
time with the length. Indeed an independent single ex- 
ponent fit of the Taute et al data [38| gives r oc L a 
with a — 2.9 in the experimental range considered 
2.2/xm < L < 28/xm. This peculiar scaling can be un- 
derstood as originating from hydrodynamic relaxation of 
a "wobbling" polymorphic chain. Using only the pre- 
viously best fitting parameters of the static data (Fig. 
11) (n /q ) 2 ?» 4. 8 x lfT 3 and rj = 10" 3 Pa • s (wa- 
ter viscosity) we find a remarkable correspondence be- 
tween the theoretical prediction Eq. [32j i.e., r t h (L) « 
L£ ext {L)/k B T and data of [H], as shown in Fig. 12. This 
0-parameter prediction matches well the data for larger 
lengths. For scaling comparison, it is interesting to com- 
pare the data with the approximate theoretical relaxation 



time Eq. 



r(L) 



3k B T 



{ k o/qo) L (expression 



strictly valid for L ^> 0.8 /im) which has a scaling law in 
agreement with the single exponent fit of the data. Con- 
sidering that £j_ ~ I.677 — 2.3?7 is roughly length indepen- 
dent in the experimental L range one can assume £j_ ~ 277 
and obtains the prefactor T t h/L 3 — 7.9 x 10 14 s/m 3 . Keep- 
ing in mind the simplicity of the interpretation (and 
the lack of free parameters therein), this compares very 
favorably with the best fit of experimental data slope 
T ftt /L 3 = 6.25 x 10 14 s/m 3 (cf. Fig. 12). Note that 
the approximate value of Tth seems to correspond a little 
better to the data but this is likely of no physical sig- 
nificance at this level of approximation. Indeed the ne- 
glected elastic modes other than the wobbling as well as 
the approximation of the hydrodynamic friction should 
change the details of the L dependence of £ ex t and thus 
of T t h (but not the general trends). Without a more 
precise computation of the dynamics we can be satisfied 
with the rather astonishing agreement between theory 
and experiment for long MTs. This leads us again to the 
conclusion that in these experiments long MTs behave as 
almost rigid helical polymorphic rotors whose motion is 
dominated by the zero energy "wobbling" mode and its 
hydrodynamic dissipation [62j. 

For very short MTs we should expect deviation from 
this simple interpretation. In this regime the linearly 
scaling internal dissipation, coming from the migration 
of polymorphic defects, should start to dominate over 
pure hydrodynamic friction and for sufficiently short MT 
lengths L — > we could measure £j„t from the limit value 
of Tth/L. It appears that for the presently available data 
L > 2 /xm [38j | this plateau-regime is not yet fully de- 
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veloped, enabling us to provide only an upper numerical 
estimate for the inner dissipation £j„t <4x W~ 17 Ns. 



comparison of our theory with the passive gliding assay 
appears more appropriate and gives a closer agreement. 
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FIG. 12: The experimental microtubule relaxation times [38| 
and the no-adjustable-parameter theoretical prediction (full 
line) as obtained from static data in Fig.ll (with Ib = 25mm, 



X — 7.5/im, 



18/im, qol^ 3> 1 ). The dashed line illus- 



trates the long length approximation Eq. 
characteristic cubic scaling with length. 



Very Short MTs 



1331 displaying the 



Comparison with experiments for even shorter chains 
(L < 7r/qo) is more difficult due to the lack of data in 
clamped MT experiments (L > 2/j.m). We can never- 
theless try a comparison with the results by the Dckker 
group for the kinesin motor gliding assay of short MTs 
3J, |35| . Besides some similarities with the clamped fluc- 
tuating MTs, there are a number of differences between 
the modelled situation of free 3D MTs and the 2D glid- 
ing assay. In particular the 2D geometry will strongly 
perturb the preferentially 3-dimensional helical ground 
state. The active contribution of strong motor forces on 
the trajectory of short MTs is an additional potential 
perturbation. Effects of MT buckling and axial MT ro- 
tation by kinesins become likely important. This said 
and ignoring the differences we can still compare (in or- 
der of magnitude) our and Dekker groups' results 34 , 3|| . 



For micron sized MTs we obtain l p ~ 0.8/j,m in approx- 
imate agreement with ~ 0.2/im obtained in | 3 ^ . One 

give 



should mention that the two cited studies 
strongly different results depending on whether free glid- 
ing or gliding assay with additional electric field were 
considered. This difference comes from the larger de- 
formations induced by the electric field. Therefore the 



SUMMARY 

In summary, we have suggested a new model that con- 
nects some of the most persistent and confusing experi- 
mental findings concerning microtubules. Starting from 
a rather broad spectrum of (apparently) disconnected 
observations we have progressively built the case for a 
new hypothesis : the existence of an internal switching 
of the GDP tubulin dimer within the microtubule lat- 
tice. Why do microtubules become helically wavy, why 
do they switch to permanently bent circular states, why 
do they fluctuate anomalously when clamped? These 
three dangling questions became the central pillars for 
the present model. Surprisingly, the simple assumption 
of a bistable GDP-tubulin seems to explain these other- 
wise disparate phenomena in a unified manner. As we 
know from recent experiments - the bistability hypothe- 
sis of taxol stabilized protofilaments is an empirical fact 
indeed [25[. We have shown here, that the incorpora- 
tion of such a bistable tubulin into a closed elastic lattice 
changes its free behavior - it introduces strong conforma- 
tional competition among the tubulin dimers. Tubulin 
units on opposite sides of the tube now start to compete 
for who is going to switch to the curved state. The lattice 
induced frustration does not allow all the tubulin dimers 
to minimize their energy individually and to switch to 
their preferred states at the same time. The symmetry 
breaking induced by this frustration mechanism leads to 
a global microtubule lattice curving. Remarkably the 
curving direction is chosen randomly - and this has pro- 
found consequences. The microtubule can chose between 
many energetic ground states (as many states as protofil- 
aments in the lattice) . When we graft one end of the mi- 
crotubule onto a substrate while still allowing it to chose 
its bending direction freely the strange energy degener- 
acy generates a very unusual thermal motion. In this 
case the microtubule's motion follows -roughly speaking- 
a cone and it rotates - or "wobbles"- at no energy cost 
around its attachment tangent. This mode of motion - 
which is not to be confused with material frame rotation 
(which is strictly prohibited by grafting) - is probably 
among the most striking outcomes of the two state GDP- 
tubulin model. It is exactly this behavior that allows to 
consistently explain the measurements of unusual lateral 
fluctuations of grafted microtubules (37 , 



PERSPECTIVES 

We have focussed here on modelling taxol stabilized 
microtubules and the question naturally arises if the 
model proposed here affects the 'real' in vivo micro- 
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tubules. On the one hand the 'weakly curved' state 
that is involved in the soft polymorphic dynamics as de- 
scribed here seems to be (so far) a specific signature of 
taxol stabilized GDP tubulin state. On the other hand 
we have argued that the naturally occurring 'high curva- 
ture' GDP-tubulin state could coexist with the straight 
state in the lattice under in vivo conditions where MTs 
are stabilized with MAPs. The involvement of this 'high 
curvature' state switching seems to manifest itself in mo- 
tor driven straight to wavy transitions of MTs in many 
living cell systems [3^, Eil - tH^ . A particularly impres- 
sive instance of such a polymorphic switching event in 
vivo could be found in the process of axonal retraction 
where the whole MT cytoskeleton of the axon undergoes 
a straight to helical transition and in turn retracts to- 
wards the soma [Hij . To understand these dramatic tran- 
sitions in vivo the present theory has to be advanced and 
modified in two manners. First, the effect of large ac- 
tive motor forces (rather than thermal ones) has to be 
taken into account. In particular, one expects that un- 
der strong buckling forces even thermodynamically unfa- 
vorable states can become activated and constitute the 
ground state upon large loads. 

Second, in virtually all in vivo experiments MTs are es- 
sentially confined in 2D as the containing cells adsorb to 
the glass substrate and assume a vary flat 'fried egg' con- 
figuration. Consequently, the measured properties will 
not necessarily reflect the three-dimensional properties 
of the molecule. This is particularly important for a MT 
transformed to a polymorphic helix state where the con- 
finement entails naturally a strong deformation (of the 
initially 3 dimensional ground state). Under confinement 
the helical bending and torsional modes become strongly 
coupled and bring about new physical effects. In partic- 
ular a torsionally very soft helix will have a tendency to 
unwind and form in extreme cases circular arcs - remi- 
niscent of the rings observed in gliding assay experiments 
42]. 

Finally, the local action of molecular motors could trig- 
ger a switching to the highly curved state. While for 
classical motors like kinesin 1 direct evidence for such 
a mode of action is still missing its relative kinesin 13 
[25L l58| has a well documented ability to actively trig- 
ger radial bending of protofilaments. Other molecules 



like katanin 

Hi |60( have also been suggested to perturb 
the lattice and trigger longer range transitions 46] . This 
opens the intriguing question : could classical motors (ki- 
nesin 1 and dynein) trigger cooperative state transitions 
and even transmit conformational signals along the tube? 
Considering the present model for stabilized MTs (where 
high cooperativity is inherent to the data interpretation) 
this idea might not be far fetched. In fact some evidence 
towards long range cooperativity of kinesin binding along 
the MT was presented by Muto et al. |6l| - however these 
results still await robust reproduction. 

This brings us to the question of what experiments 



should be performed in order to nail down the polymor- 
phic mechanism or any other mechanism for MT dynam- 
ics. With microtubules being such delicate, subtle and 
possibly long range correlated objects (as suggested here) 
a general rule of thumb for experiments should be: Treat 
them more gently (do not confine) and observe more care- 
fully (look for correlated motion). A simple yet impor- 
tant experiment would be a systematic direct observa- 
tion of one-side grafted but otherwise completely uncon- 
fined MTs fluorescently labelled along their full contour 
length. As mentioned the presence of a quasi 2D confine- 
ment in thin chambers as used in most MT experiments 
so far would perturb the native helical MT state and 
should be therefore explicitly avoided. The freely sus- 
pended gold-EM nanogrid attachment geometry as used 
by Pampaloni et al. seems particularly suited for that 
task. Going beyond Pampaloni et al. who labelled and 
traced the MT end only (via a bead), the microtubules 
should be visualized along the full contour in this ge- 
ometry. Tracing of several or all points along the con- 
tour should reveal the predicted sinusoid- helical nature 
of MT states. The present model predicts a peculiar 
cooperatively rearranging helix state with characteristic 
tell tale curvature correlations between different lattice 
positions which are entirely absent for usual semiflexible 
filaments. Directly observing such collective motions - 
like the suggested "wobbling mode" - while prohibiting 
trivial spacial rotations that could mask the effects (by 
MT grafting) would constitute smoking gun evidence for 
a polymorphism related mechanism. 

In conclusion, we have proposed a novel model for 
internal MT lattice dynamics. We have shown that it 
accounts for the otherwise mysterious MT helicity [27j |. 
the anomalous length dependent lateral fluctuation static 
[37|,|38| and dynamic scaling [38J. The latter two phenom- 
ena appear as mere consequences of the peculiar "wob- 
bling motion" of the polymorphic cooperatively switching 
MT lattice. Although most of the observations discussed 
here are made in vitro on taxol-stabilized MTs, we pro- 
vided arguments in favor the existence of polymorphic 
MT states in vivo. We speculate that the implied con- 
formational bistability of tubulin and the allosteric in- 
teraction are more than just nature's way to modulate 
the elastic properties of its most important cytoskeletal 
mechano-element. It could also be a missing piece in 
the puzzle of polymerization " catastrophes" . Even more 
intriguingly the predicted structural cooperativity could 
allow for long range conformational signalling along sin- 
gle MTs and turn the latter into an efficient " confotronic" 
wire transmitting regulatory signals across the cell. 
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APPENDIX 
A. Polymorphic phase coherence length 

In this section, we derive the formula — 
-g^y (2 + e 2J / k BT^ f or £ ne p iy mor phi c phase coherence 
length. To this end we want to calculate the distribu- 
tion of double junctions that leads to angular orientation 
change A<I> on a scales I much larger than the tubulin 
dimer 6, yet still much smaller than the total length: 
b <C I <C L. In this domain, at each cross section we 
have 3 possibilities: 

1) State j = with no double defect. The rotation 
angle A$ is attached to the internal lattice rotation, 
qo =0 



A* 
b 



2) State j = —1 for a left handed double defect, A$ 
deviates from the internal twist : A£ _ q Q = _ I 2jl 

3) State j = +1 for a right handed double defect with 

A<S> „ _ , 1 2tt 

On a length I we are throwing a 3 sided dice l/b times 
and the total rotation of A$ away from optimal twist is 
A$— q a l = ^ J2nti in- The variation of the polymorphic 
phase with respect to the internal twist is then 



A$ 1 2tt 

A0= — -90 = 7-^2^" 



71=1 



For I /a » 1 the law of large numbers implies that the 
random variable A<j) = jj^J2n=iJn becomes Gaussian 



distributed 



p (A<j>) oc e 2 < A * 2 > 



with mean (A</>) = and (Acj) 2 ) = (£) (j 2 ) (as 

(jnjm) = <$mn(j 2 ))- The average (j 2 ) is given from 
the Boltzmann factors of the three different states po = 

l+2e-^ J a ^ ±:L = l+2e- 2 ' llJ SO that (j ) = l+2e -ifJJ ■ 

We can now interpret the quantity 1/(2 (A0 2 )) as com- 
ing from an effective elastic energy over the interval I by 
writing 2 {A<t> i ) — \PC ( ) > 1 (A<j>) 2 which allows us to iden- 
tify the effective stiffness 



8^ 2 



C 6 = fcT^ (2 + e 2Jfj ) b. 



Note that this expression is valid for large enough 
J suppressing higher order defects i.e. in the limit when 
multiple double defects sitting on a single lattice site (i.e. 
\j\ > 1) can be ignored. 



B. The variation of the polymorphic modulus 

In this appendix we compute the energy variation due 
to a deviation of the polymorphic modulus \P\ away from 
its optimal value \P*\ minimizing the energy; i.e., the 
change of the number of switched PFs. We start with 
the energy density of a MT cross section 



e = f K P oi{p)f + n 2 (7-^5 -sin 2 (jjp))) 

(34) 

whose minimum energy is reached for p* = y — 
A- arcsin7. Assuming a continuous number of PFs, the 
energy of a state with p = p* + Ap switched PFs reads 
to quadratic order: 



e{p* + Ap) « e(p*) - cos C^p* 



Ap 2 



e{p* 



ir 2 B 
N 2 



kW1-j 2 A p 2 



where we used cos(7r — arcsin7) = — yl — 7 2 - Therefore 
the energy variation of a segment of length I reads 



AE 



n 2 B 
N 2 



i4^1^ff ds{p{s)-p*{s)f 



Now using |P(s)| = |sin(^p)| / sin(n/N) we can write 
the energy variation to the same (quadratic) order as 

AE « BKls\n 2 {-K/N)jl^ff Q ds{{\P{s)\ \P*\)f 

Therefore any deviation of P from its optimum state 
P* is associated with an energy cost proportional to the 
length I of the region in the unfavorable state. 



C. Persistence length(s) 

A definition of the persistence length, often used 
in single molecule experiments, is expressed in terms 
of the lateral deviation ~f? = (x(s),y(s)) of a MT 
clamped at s = from its attachment axis : I* (s) = 

2/3s 3 /((7 > ( s ) - ("^( s ))) 2 ) and (••> is thc statistical 
average. The equivalence of thc x and y directions 
implies that l* p (s) = l/3s 3 / (t(y (s) - (y (s)» 2 ^ . Thc 

second often used alternative but more standard def- 
inition of the persistence length - the tangent per- 
sistence length- is related to the angular correlation 
l p (s — s') = \s — s'\ /V(s — s') with the variance V — 

(i&v ( s ) — ®y ( s ')) 2 ^ (by symmetry wc have the same ex- 
pression with 6 X ). Whereas for an ideal WLC I* — l p = 
Ib is position and definition independent this is not the 
case for a polymorphic chain (see Fig. 13). For small an- 
gular deformations the decoupling of chain's fluctuations 
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into polymorphic and purely elastic contributions allows 
to decompose the persistence length as l~ x = + 1^ 1 - 
this result being valid for both definitions of the persis- 
tence length. 

We first focus on the first definition - the clamped per- 
sistence length. In this case the polymorphic persistence 
length lt i(s) is given by 



(35) 



i; ol (s) = l/3s i /((y pol (s)-(y pol (s)))' 



where y po i(s) is the lateral polymorphic displacement in 
the y direction. Integrating over the rotational zero mode 
readily implies (y pol (s)) — (see Eq. [TS]). From Eq. [Tj5] 
one can write 



{yloii. s ))=l \ G(si, s 2 )dsids 2 
Jo Jo 

with the angular correlation function G(si,S2) — 
(9y,poi{si)0 y , P oi(s2)) given by the integration over the 
zero mode 



(36) 



of the angular correlation function at fixed value 
of 0o, i-e., G (si,s 2 ,4>o) = {dy, P oi(si)9 VtPO i(s 2 )) \<j, a . 
This last expression, from the relation 9 VtPO i{s) — 



k f* sin (4>(s') 
itly given by 



qos 



o ) ds' (cf. Eq. [HI) is explic- 



Go(si,s 2 ,(f>o) = kI J J (sm((t>(s) + qos + <po' 

sin (j) (s') + q s' + </> ) ^ \<f, dsds'. (37) 

After integration over 0o and using the known re- 
sult {^cos (Z(si) - 0(s 2 ) n = e -l s i- s 2l/ 2 '* which results 

from the WLC type probability distribution of the field (f>, 
i.e., P[<f>] ~ exp(— l -j Jq dscj)' 2 ) one obtains the rotational 
invariant correlation function in the form 



2 r s i r s % 



G(si,s 2 ) = ^/ / e~ 21 « cos (q (s-s'))dsds' 
2 Jo Jo 

(38) 

Computation of the integrals in Eq. [35] gives finally the 
following expression for the polymorphic contribution of 
the transverse displacement 



(v P oi{s) 2 ) 



+e 2iS *F 3 sin(g oS )} (39) 
with Pi(s)= 24Z| (l - 6x + x 2 ) -3Z (l + x - x 2 - a; 3 ) s 2 
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FIG. 13: Different definitions of the persistence length can de- 
viate from each other for a polymorphic chain. The "clamped 
persistence length" lp (thick line) and the "tangent persis- 
tence length" l p (thin line) (for Ib = 10mm, 1$ = 50/im, 
fto = 0.03/im -1 and qo = 0.7/im -1 ). 



+ (l + 3x + 3a; 2 + x 3 ) s 3 , P 2 (s)= 24Z| (l - 6x + x 2 ) 
+I2ll (1 - 2x - 3x 2 ) s and P 3 (s)= 192/^(1 - x) 

+24/^(70(3 + 2x — x 2 )s where we have introduced the 
notation x = il^q 2 . 

From Eq. 1391 we get the polymorphic persistence 
length l* ol (s) defined in Eq. [351 and in turn the global 
persistence length l*(s) depicted in Fig. 13. Its physical 
interpretation is discussed in the main text. 

We now consider the second definition of the persis- 
tence length l p (s - s') = \s - s'\ /V(s - s'). From Eq. 
1381 the angular variance V po i can easily be evaluated 



V po i(s) 
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>(i-4gg$)' 

1 + tilql 



+ 



^((l-4^J)cosM 
-4(7o^0 sin (q s)) 



(40) 



The resulting persistence length l p (depicted in Fig. 13) 
shows a rich behavior similar to the persistence length 
I* (s) but displays a distinct functional form from the 
latter. However as expected, both curves reach the same 
asymptotic value at very short and very long MT lengths. 



D. Zero mode dynamics 

The evolution of the zero mode <j>o (t) is given by Eq. 
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al t,tot JO 



(41) 



with a friction constant f tot = f int + £ e:Ef where £ ext is 
given by Eq. 1271 The correlation function of the ther- 
mal white (s, t) noise is (T^ (s, t) (s 1 , t')) — DS(s — 
s')8(t — t') with a coefficient D that can be determined 
in the following manner. Notice first that performs a 
free Brownian motion and its quadratic fluctuations nec- 



essarily satisfy the relation \ {4>a (t) — 4>q (0))' 
On another hand integrating Eq. [U : 



2k B T 



00 W ~ ^o (0) = S 



*jf (zjf r(a, *' )da ) <ft ' 



(42) 



and averaging over the white noise the quadratic phase 



fluctuations one obtains 



(0o (t) - 0o (0))' 



- D / 

from which we readily deduce I? = 2^ tot fcsT - as ex- 
pected from the fluctuation dissipation theorem. 

The relaxation time is generally given from the time 
correlation function < y P oi(s,0)y po i(s,i) > with the 

lateral position y po i(s,t) = ^ (sq cos (<j> (t) + a) + 



sin (0o (t) + a) — sin (qos + 0o (t) + a)) obtained from Eq. 
[14] with » s. The average must first take into ac- 
count all statistically equivalent values of angular ori- 
entations a e [0,27r], such that (y po i(s,0)y po i(s,t)) — 



Jo* {y P oi(s,0)y po i(s,i)) a |f and we obt 



am 



(y P oi(s,0)y po i(s,t)) = (yl i(s)) (cos(0 o (t) - 0o (0))} 



with<^ (s) >= |(V + I: 



■cos(ijos) 



(43) 
^ LJ - cor- 

90 / 

1 in the limit 



responding to the static result Eq. 
s/lij) « 1. With W2\ defining a simple Gaussian random 
walk processes one straightforwardly obtains 



(cos(0 o (t) - 0o (0))) = e-*/ T 
with the relaxation time given by 

Ctot 



t = L- 



kpiT ' 



(44) 



(45) 



E. Comment on MT surface attachment and the 
robustness of "wobbling" 

Throughout this work we have assumed that the free 
rearrangement of the polymorphic lattice states is not sig- 
nificantly hindered by the covalent surface attachment of 
the MT, as e.g. performed in (37| and (3^ . This assump- 
tion is integral to the "wobbling" motion and in turn to 
understanding the static and dynamic data scaling. It 
therefore deserves a closer consideration. 



In the experiments by Pampaloni 37[ et al. and Taute 
et al. [38| the adsorbed MT part is attached to a gold 
(electron microscopy grid) surface via thiol groups. It is 
likely that ks 1 - 2 protofilamcnts will establish localized 
chemical contacts with the gold microgrid. While a sub- 
stantial perturbation of the dimer like e.g. denaturation 
appears unlikely, it is unclear to what extent this pro- 
cedure will perturb the inner (polymorphic) dynamics of 
the entire tubulin dimer units. In principle one can an- 
ticipate two plausible scenarios that would to a varying 
degree interfere with the free "wobbling" motion: 

51) Due to high cooperativity (large coupling J) the 
polymorphic state transition can propagate within a 
certain penetration depth into the adsorbed (straight- 
planar) MT section. 

52) The cooperativity is too weak to compete with the 
constraints imposed by the surface (including chemical 
perturbations) and the polymorphic transition does not 
propagate into the straight adsorbed MT section. 

In both cases we have a non vanishing deflection an- 
gle between a forced (adsorbed) planar section and the 
free helical section direction - causing effectively the char- 
acteristic MT "kink" at the surface interface. However 
the rotational mobility of this "kink" (wobbling mode) 
which is integral to our theory will be affected in slightly 
different manner. 

If in case SI in the adsorbed section the polymor- 
phic order parameter P can rearrange to some extent 
(by switching the monomer states without causing a de- 
tectable deformation) except for possibly in the few sur- 
face interacting dimers, then the effects of the "wob- 
bling" motion will be hindered only mildly in the follow- 
ing sense. To retrieve the anomalous lateral fluctuations 
it is indeed enough for the wobbling angle 0o to move 
freely in a certain non- vanishing angular interval. A sin- 
gle complete or multiple rotations of the order parameter 
P are not strictly necessary for the "hinge" effect - and 
they are in fact equivalent in lateral projection (as in ex- 
periment) to the motion of the wobbling angle 0o in the 
smaller interval [— tt/2, +tt/2]. Note that even smaller 
intervals than that will lead to a similar phenomenol- 
ogy (in particular dynamic and static variable scalings 
with length). Thus the conical hinge- like motion is in 
a sense robust with respect to a limited local rotational 
hindrance perturbation in the adsorbed region. 

In the scenario S2 the situation is somehow simpler as 
the polymorphic dynamics of the adsorbed region is not 
involved in the process (the polymorphic order parameter 
vanishes there: P = 0). Wobbling is realized through a 
coherent rearrangement of the free MT section alone - 
without a strong coupling to the adsorbed region. 

Although both attachment scenarios SI and S2 appear 
to some extent plausible, at present it is difficult to make 
reliable statements about their respective likelihood. In 
fact only a posteriori we can cautiously state that based 
on the experimental static and dynamic measurement ev- 
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idence: the chain "wobbles" to a high enough extent to 
display the effects that we observe. 
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